Generation of field mediated three qubit entangled state shared by Alice and Bob 
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A scheme to generate shared tripartite entangled states, with two-trapped atoms in a cavity 
held by Alice (qubits Ai and A2) entangled to a single trapped atom in a remote lab owned by 
Bob [B), is proposed. The entanglement is generated through interaction of trapped atoms with 
two mode squeezed light shared by the two cavities. The proposed scheme is an extension of 
the proposal of ref. [W. Son, M. S. Kim, J. Lee, and D. Ahn, J. Mod. Opt. 49, 1739 (2002)], 
where the possibility of entangling two remote qubits using a bipartite continuous variable state was 
examined. While the global negativity detects the free entanglement of the three atom mixed state, 
0^ ' the bound entanglement is detected by the negativity calculated from pure state decomposition 

^^ ' of the state operator. The partial negativities calculated by selective partial transposition of the 

^^ , three atom mixed state detect the pairwise entanglement of qubit pairs A\B, A2B, and AiA2. 

^^ ■ The entanglement of three atoms is found to be W-like, no GHZ like quantum correlations being 

r^ ' generated. 
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INTRODUCTION 



f^ . Quantum systems with discrete energy spectrum are often used for storing, manipulating, and transmitting infor- 

I ' mation encoded in quantum states belonging to finite-dimensional Hilbert spaces. A two-level atom is an example 

r-| I of two-dimensional quantum systems, universally known as qubits. All possible superposition states of a two-level 

C^ ' system can be expressed in terms of linear combinations of two-dimensional orthonormal states defining the standard 

^ I computational basis. On the other hand, the electromagnetic field is characterized by continuous variables that is ob- 

^^' servables with a continuous spectrum of eigenvalues. Quantum information based on continuous variables is encoded 

in quantum states belonging to infinite-dimensional Hilbert spaces. The entanglement of remote qubits is a resource 

^s) . for implementation of communication protocols [1] and information processing [2] using quantum systems. W. Son 

^ ' et al. [3[ have shown the possibility of entangling two remote qubits using a bipartite continuous variable state. In 

^^ . this article, a scheme to generate shared tripartite entangled states of two-trapped atoms in a cavity held by Alice 

2l ' ( qubits Al and A2) and a single trapped atom in a remote lab owned by Bob (B), is proposed. A single mode of 

^^ [ the two-mode squeezed electromagnetic field, generated by a squeezed light source, is injected into each cavity by a 

• ■ beam splitter. The entanglement is generated through interaction of trapped atoms with quantized electromagnetic 

f^ ' field shared by the two cavities. No direct interaction amongst the atoms takes place. 

The dynamics of three atom entanglement due to atom field interaction is investigated after tracing over the 
field degrees of freedom. We use the global negativity and partial K—way negativities to analyze the entanglement 
dynamics. Negativity [4|, based on Peres Horodecki |5J, l6| criterion, has been shown to be an entanglement monotone 
[31 . A comparison of global negativity of three atom mixed state with respect to qubit B and the linear entropy 
calculated from the state operator for qubit B indicates that for certain ranges of interaction parameter values the 
5^ global negativity fails to detect the entanglement amongst the atomic qubits. Using pure state decomposition (PSD) 

Cu [ of mixed three atom state, we recalculate the negativity by making a weighted sum of pure states contributions 
to global negativity. The global negativity calculated from pure state decomposition (PSDG) of the state operator 
is not a proper measure of entanglement because it is not invariant with respect to local operations and classical 
communication (LOCC). It does, however, detect the bound entanglement being greater or equal to the convex 
roof extension (CRE) of negativity shown to be an entanglement monotone ^Sj. We use the global negativity and 
partial K—way negativities of pure state decomposition, primarily, to understand how quantum coherences present 
in the three atom state, evolve with interaction time. The coherences of a multipartite composite system having N 
subsystems can be quantified by partial K—way {2 < K < N) negativities @, [ifll- Partial i^— way negativity is the 
contribution of a i^— way partial transpose to global negativity. In a three qubit state, the entanglement arising due 
to 3— way and 2— way coherences |ll| is quantified by partial 3— way and 2— way negativities. For canonical states, 
positive valued partial 3— way and 2— way negativities measure, respectively, the three qubit GHZ-like and pairwise 
entanglement of the system [ij] . We further use the partial negativities calculated from the pure state decomposition 
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of three atom mixed state to study the dynamics of two-body coherences of qubit pairs AiB, A2B, and ^1^2- The 
qubits Ai and A2 are shown to form a beU pair with certain probabihty, which is entangled to qubit B. The degree of 
entanglement depends on the atom field coupling strength, squeeze parameter, as well as the coefficient of reflection 
of the beam splitter. 

II. THE CAVITY FIELD 

Consider two cavities prepared in the vacuum state at i = 0. The beam splitter used to inject the external field 
into a cavity is represented by the operator 



B{9) = exp 



^(cp-c^f) 



(1) 



where c ( / ) and c' [Pj are creation and annihilation operators for photons in the cavity (external field) , respectively. 

The coefficient r = cos | is the reflection coefficient of the beam splitter. The two mode squeeze operator acting on 
two mode vacuum state 

S{s) |0, 0) = cxp(-sa6 + sa^P) |0, 0) , (2) 

produces a two-mode squeezed field in composite state 

-. 00 
|*F> = ^^^^(tanhs)"|n,n>, (3) 

n=0 

where s is the squeeze parameter. The bosonic creation and annihilation operators for mode one are a\ a and for 
mode two are b^ , b, respectively. The two mode squeezed state is an entangled state having bipartite entanglement 
determined by the value of squeeze parameter s. Using beam splitter (represented by operator of Eq. ([T])) to inject 
the modes one and two of the squeezed field (Eq. ([3])) into independent cavities ci and C2, respectively, the cavity 
field state at i = reads as 

, ^ 2 00 min[n,m] 

where 



V cosh s , 

^ ^ n,m—Q k,l—0 



airm = c? (») C7 (0) c," (9) cr (») . (5) 



='"'*' = m^'.m^r'''- (i j ^"'""" li j "^^' U 1 ■■"^""' I ^ I ■ '^' 



and 



^'"<''> = \/s5rr«T»^H2J^'-""l2l- <" 



The operator PciC2{^) represents a mixed composite state of cavity field. 



III. ATOM FIELD INTERACTION 



The interaction of N two-level atoms with a resonant single mode quantized electromagnetic field in rotating wave 
approximation is described by Tavis -Cummings model [l3| Hamiltonian 

i=l,N i=l,N 




m: 
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FIG. 1: The schematic diagram of experimental setup 



where luq is the atomic transition frequency and g the atom field coupling strength. Eigen states of Pauli operator ct* 
model the internal states of the i*'* atom {i — 1, 2) with eigenvalue m^ = —1 (+1) standing for the ground (excited) 
state of the atom. Defining the composite operator, at — '^iO'l. where k = (z,+,— ), we may construct the eigen 
basis of operators ct^ and ct^ to represent N atom internal states. A typical coupled basis vector is written as |cr, Too-), 
with eigenvalues of ct^ and ct^ given by a{a + 2) and TTig-, respectively. 

Consider two atoms trapped in cavity ci radiated by the single mode resonant cavity field in photon number state 
|n). Interaction with the cavity field generates entanglement of internal states of the atom, and the cavity field state. 
Working in the atom-field basis \a, m^,, n), the relevant basis states for the case when two atoms are prepared in their 
ground states at t = are |2, —2, n), |2, 0, n — 1), and |2, 2, n — 2) and the Hamiltonian matrix reads as 



hgBn 
H^i = I figBn hgA, 
hgAn 



(9) 



where Ar. 



exp 



v/2(n-l), B„ 



-[iHut\/n 



Ui{r) 



2n. Defining interaction parameter r — gt, the unitary operator Ui (r) 



is found to be 



^ [Bl cos (/„r) + Al] -iB„ sin (/„r) A„B„ [cos (/„r) - 1] N^ 

Al + Bl ^(Al + Bl) Al + Bl 

-JBn sin (/ut) , , —iA„ sin {fnr) 



^{Ai + Bl) 

A„B„ [cos jfnT) - 1] 
■P^n + -On 



-lAn sin (fngr) 



ViA'i + Bl) 

[Al COS (/nr) + Bl] 

^n + B„ 



(10) 



where /„ — ^2 (2n — 1). 

For a single atom interacting with the cavity field in cavity C2 , the unitary matrix that determines the time evolution 
of the system in the basis |1,— 1,to),|1,1,to,— 1) is 



U^ir) = 



cos {y/rnr) —i sin (y^r) 
i sin (-\/tot) cos (-y/ror) 



(11) 



The evolution operator Ui2''{t) for the two cavity composite system is obtained by taking the tensor product 



t^iTM = t/rw®c^rM- 



(12) 



IV. THREE QUBIT STATE OPERATOR 

Consider two-trapped atoms in cavity ci held by Alice ( qubits Ai and A2) and a single trapped atom (qubit B) 
in a remote cavity C2 controlled by Bob. The atoms are cooled down to their ground state with cavities prepared in 
vacuum state at i = 0. A squeezed light source generates the two mode electromagnetic field. A schematic diagram of 
the setup is shown in figurc((T]). The two mode squeezed field is injected into the cavities, with each cavity receiving 
a single mode of the field. The composite state of field in two cavities is given by PciC2{0) as in Eq. fH). The state of 
the composite system without atom field interaction is represented by 



P(0) == P 



A1A2B 



(0) (g, PccjO), 



(13) 



where ^^^^^^(0) — |2, —2)^^ |1, — l)^^ (I7 ~l|ci (^7 ~2|^^ in the coupled basis (basis states labelled by eigenvalues of d^ 
and CT^). Here |2, —2)^ represents qubits Ai and A2 in state |2, —2) in cavity Ci. The state operator, after interaction 
between the atoms and the respective squeezed fields in the two cavities for a time t = r/g, is given by 



P(r) = U^2{r)p^'^^^{0)^Pc,cMUl2{r) 



(14) 



Using the unitary evolution operators due to atom- field interaction in cavities one and two (Eq. (fT2|) ). the state of 
the composite system after interaction time t is written as 



P(r) = 



/ 1 



f/l2(T) 



yj (tanhs) 



n-\-rn 



.n.m— 



\ cosh s 

n[n.rn] 

^ Gr(0)|2,-2,n-fc)|l,-l,n-O(2,-2,™-fc|(l,-l,m-/| 



mill n.m 



k,l=0 



f 1 



uUr) 



V cosh s 



2 00 niin[n,m\ 



J2 E (tanhs)"+'"G^rw{^r'w®f7r'w 



n,r?j=0 k,l=0 

|2, -2, n-k) |1, -1, n - I) (2, -2, m - k\ (1,-1, m - ;| ^/^""'^(t) (g, {/^'""'(t)} . 

The reduced state operator for atoms is obtained from p{t), by tracing over the field modes, that is 

p^i^^^(r)=Tr,,,,(p(r)). 



(15) 



(16) 



Associating, single qubit states |0) (|1)) with ground (excited) state of the two-level atom, the initial state of 
three atoms in computational basis reads as ^"^^"^^'^(O) = (|00) (00|) g) (|0) (0|). The matrix p^^^^^{t) in the 
computational basis |00),^ |0),^, |10),^ |0),^, |01),^ |0),^, |11),^ \0)[^, |00),^ |1) J |10),^ |1),^, |01),^ |1),^, |11),^ |1),^, 
reads as 
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(17) 



with the matrix elements given by 

/ -, \ 2 oo n 



cosh s , 

' n=Okd=0 

ih ~ rZ _L it / n \ \ / \ 



P2i^'''{r) = i—A 5]5](tanh.)2«G;!rWcos2(TV^^^ 



(n-k) 2 



sin' ifnkT) , (19) 



(2n - 2/c - 1 



2 oo n 

, cosh s , 

^ ^ n=Okd=0 

(n — k) (n — k — 1) , , „ ^ ,2 , x 

(2n-2fe-l)^ (-(^"--)-^)' (20) 



/ -J \ 2 00 n 

pil^^^ir) = — — ^^(tanh.)2"G^r(^)sin'(TA/^r^ 



cosh s , 

' n=Ok.l=0 



(2n-2fc-l) 



{{n^ k - 1) + (n-k) cos if nkT)f, (21) 



2 00 n 



Pti^-'^ir) = Mj^j 5:5:(tanh.)2"G;r(^)sin2(rV^r^ 



n=Ofe.(=0 



^"^ ^^ ■sin2(/„fcr), (22) 



(2n - 2/c - 1) 

/ -. \ 2 00 n 

pT'^'ir) = — j- ^ ^(tanh.)2"G;r(e)sin2 (rV^T-? 

(n — fc) (n — fc — 1) , , , , ,,2 /„,N 

(2n-2fc-l)^^ ^-^^(^"^-)-^)' (23) 



2 00 n 



pT'^'ir) - (-^) EE(t-i^^)'"^^Gir^(^) 



^o^h^/ t'okU V(2n-2fc + l) 

cos f rV?^ — ^ ) sin I rvn^^ + l) sin(/„+ifeT) 

(n— fc— 1) (n — k) ,„ ,\ ,„,, 

^ ^ ' ^ ^ cos (/„fcT) , (24) 



(2n-2fc-l) {2n-2k-l) 



,4^^^-(r) . -i- 5:5:(tanh.)-+iGir^(^) 



cosh.y ^^^^^^ V(2^-2A:"1) 



\/{n~-k){n-k+l) / ^ -\ 

— — cos TVn — I 

{2n-2k+l) \ J 

sin [r^/n - / + 1 j (sin (/„fcT)) (cos (/„+ut) - 1) , (25) 

and /„fc = V2(2n-2fc-l). 



S=1.2 




FIG. 2: The probabilities, versus r(= grjt) for s = 1.2. 



V. PROBABILITY OF GENERATING THREE QUBIT ENTANGLED STATES 

We have calculated, numerically, the matrix elements of state operator p ^^^(r) using a simple fortran program. 
The sum ^ in Eq. (j4]) is approximated by ^ , where the value of n^^x is determined by the choice of parameter 

s. By analyzing the variation of p^j^ ^ (r) (i = 1 — 8) with parameter s we have chosen s = 1.2 as the preferred value 
of squeeze parameter and rimax = 80. A direct observation of —10 dB squeezing of quantum noise of light, reported 
recently by H. Vahlbruch et al. |lj|, makes the selected s value a reasonable choice. 

Figure ^ displays the probabilities Pi (1 = 1 — 8) of finding the three atoms in the local basis states \iii2)c ^3)0 
(i = ii + 2^2 + 4^3 + 1), for s = 1.2 and 9 — tt . We notice that the probability of finding the atoms in their ground 
states drops to 0.34 at r = 0.8, which is already less than ^ P^ at r ~ 0.4. For s = 1.2, the probability of initial 

i=2-8 

atomic state interacting with the field is reasonably large, while the mean photon number is not too large. In case 
the mean photon number is large, the cavity decay must be taken into account. 
The operator p'^i'^^^(r) of Eq. pT)) . can be rewritten in diagonal form as 



2=1 

where pi, and |$i) are the eigenvalues and eigenvectors of p^^'^^^(t) given by 

1*1) = aooo|00),jO),^+aioi(^^^^^i^^)|l) 
A , ^A 
P' = - 2 



(26) 



^-^^^-y{pA^Pi.y~^Hpt.y 



(27) 



l$2) - 6000 |00),^ |0),^ + 5ioi ('^^^k:^^^ 1 ID 



V2 



P2 



{pti + pis) 



y{pA-pi,f+HPt,) 



(28) 



/IIO), +|01), \ 
$3) = cioo (^ ' '^'J '"' ) |0>,, + cm 111),, |1),, , 

— + yiP^2-pi,T+^p'2e^ (29) 



P22 

P3 = 



|$4) = dlOO (^^^^^'^^) I0)c. + ^111 Ill>c, Il>c. , 



Pi2 - "^ 
P4 = 



^ - yK2-p4)' + 4pi6, (30) 



2 2 

l'i>5) = |ll),JO),,, P5 = P^3, 

l'i'6> = |00),Jl),^, P6 = pi- (31) 

The coefficients oooo J Oioi, ?'ooOi &ioij Cioo, cm, anddioo, ^iii are known in terms of tfie matrix elements of p^^ (r) . 
Tire probabilities Pi determine the entanglement distribution amongst the three qubits in the system. The states |$i) , 
i = 1 to 4 are W-like states. The state p^^^^^(r) is a mixture of W-like states. For example, the fidelity of projecting 
out the state 

\Wi) = ^ (100),, |0),^ + |10),, |1),^ + 101),, |1),J , (32) 

is given by 

F{W,) = Tr (p^^^^^(t) \W,) {W,\) . (33) 

Alice can perform a controlled not gate on qubits Ai and A2 transforming the state into a mixture of GHZ like states. 

It is easily verified that the probability of projecting out the qubits A1A2 in state ( ^^^ — — j from the reduced 

state p^i^=(r) = Tre (p^i^^^(r)) is {pti + pis) ■ 

VI. THREE QUBIT MIXED STATE ENTANGLEMENT DYNAMICS 

The global partial transpose [4] of p^^'-^^^lr) with respect to sub-systems, Ai, A2, or B is obtained by using the 
relations 

{^li2^3\ P^g' (r) busjs) - (ji*2*3| P''^-'^^(t) |zij2j3) , 

(n*2*3| pI'^ (t) bU2j3) = (hJ2*3| P-'^-^^^(t) b-l*2j3) , 

(^l^2^3 1 Pg" (^) \h]2]z) = (Ji*2j3 1 p-^^^^^ (r) \nj2iz) . (34) 

The i^— way partial transpose [ij] {K = 2,3) of tripartite state p-^i-^^s t^--^^-\^ respect to qubits Ai, ^2, or B is 
constructed by applying the following constraints: 

3 
{n^2^z\ /k' IJ1J2J3) = {n^2^^\ P^^^^^ N1J2J3) , if E (^ " ^»">.^™) = ^' 

rn— 1 
3 

(n*2«3| P^"^ IJ1J2J3) = (*1*2*3| P^^^^'' IJ1J2J3) if E (1 - '^^"..^".) ^ ^' (35) 



m— 1 



(n^2^3| P^'=^ IJ1J2J3) = (*1J2*3| P^^^^^ |J1«2J3) , if E (1 - ^»".^^") = ^' 

m— 1 
3 

(n^2^3| P^"^ IJ1J2J3) = (*i*2*3| P^^^^^ b-ij2j3) if ^ (1 - 5,^,jJ ^ K, (36) 



and 

3 

{iii2iz\ p]f \3inn) = (*1*2J3| P-^^-^^"" b-ij2«3) , if ^ (1 - 5,^,, J = K, 

3 

(^l^2^3| pI" IJ1J2J3) = (*i*2*3| P^'^'"" IJ1J2J3) if ^ (1 - S,^,,J ^ K, (37) 

rn— 1 

where Si^^j^ = 1 for i^, = jm., and (5^,^^J^ = for i,^ ^ jm- 

The negative part of p^ {p = Ai or A2 or i?) contains information about the amount and nature of entanglement 
present in the tripartite system. By using a ii'— way partial transpose one can extract the quantity of K—way 
coherence (that is the quantum correlations present in a A'— qubit GHZ like state) in the three qubit state. The 
contribution of a X— way partial transpose to the global negativity is calculated by noting that the global partial 
transpose with respect to qubit p, may be written as 

pS" = ^^ + P2' - P- (38) 

pQ ) = 1 1 the negativity of Pq is given by 

iVg = -2Y,{^>f-\fj \^r) = -2E^^' (39) 

i i 

where Af "and |^f ~) are, respectively, the z*'' negative eigenvalue and eigenvector of pj . Substituting Eq. ([38)l in 
Eq. dMD, we get 

N 

-2E^r = -2EE(*M?^1*r) + 2E(*M?l*^)- (40) 

i K=2 i i 

Defining partial AT— way negativity E^ [K = 2 to N) as 

£;?, = -2E(*f1p5'|*r), while i?„^ = -2E(*f1p|*r), (41) 

i i 

we may spht the global negativity for qubit p as 

N^^^El + E^-El (42) 

The necessary condition for an A^— partite pure state not to have iV— partite entanglement is that the global negativity 
is zero for at least one of the subsystems that is Nq = 0, where p is one part in a bipartite split of the composite 
system. 

The global negativity of partial transpose of ^^^"^^^(t) with respect to B is found to be 

nB = -2 (Ar + A^-) , (43) 

where 
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and zero otherwise. Similarly 
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FIG. 3: The global negativity Nq , linear entropy 5*;, and negativity Npgj^Q, versus r(= grit) for 



= 1.2. 
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and Ao 






0, if 
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+pf/^^= 



< 



> 



„SAiA2 
l-'22 



P44 



4(pf/^^= 



/( 



„SAiA2 
P22 
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P44 



+4(pf5^^^=) 



(45) 



(46) 



Furthermore, E^ — 0, implying that the composite system state does not have three qubit GHZ state like coherences. 
The reduced state operator for qubit B is given by 

r{r)^TrA,AAP^'^'''ir))- 
The purity of a state in Hilbert space of dimension d is measured by the linear entropy defined as 

d 



(47) 



Slip) 



1 



{I - trip')) 



(48) 



The linear entropy is zero for a pure state and one for a maximally mixed state. Calculated SP Cp^iT)) and Nq are 
shown as a function of parameter r i=gt) in Figure ([3]) for s = 1.2 and 9 = tt. Whereas mixedness of the state P^(t) 
depends on the entanglement of qubit B with atomic qubits and cavity field, the negativity of p^f (r) is a measure 
of entanglement of qubit B with qubit pair A1A2. We notice that for certain ranges of parameter r, Pq' is found 
to be positive while SP value indicates that the state is entangled. Can we conclude that for Nq = 0, the qubit B 
is entangled only to the cavity field? Since the global negativity is known to detect only free entanglement of mixed 
states, the answer is no. To detect the bound entanglement, we use the pure state decomposition iPSD) of p^^^^^^ir) 
given by Eq. (P5|) . and calculate 



Nl 



PSDG 



ip)^Y.P^Nlip,)., 



(49) 



where p stands for one of the three qubits and subscript PSDG refers to global negativity calculated for a pure 
state decomposition. In ref. [g] the convex roof extension of global negativity iCRE) has been shown to be en 



10 

entanglement monotone capable of detecting bound entanglement. The convex roof extension of negativity for a 
mixed state is defined as 

N^RE ip) = min f ^ P'^G (a) ) . (50) 

where Pi = |\l/i(r)) (^,j(t)|, |^j) being the eigenvector of p corresponding to eigenvalue p^. The negativity A^p^^q (p) , 
satisfying the relation 

N^PSDGiP)>N'cREiP)^ (51) 

is not an entanglement monotone, but serves to detect bound entanglement. Numerically calculated Npgj^Q for the 
state p^^^^^{t) is also plotted (dotted line) in Figure ([3]). A comparison with Nq plot shows that the regions in 
which Nq = 0, NpgjjQ is found to be nonzero, implying that it is possible to project out entangled states from 
-^AiA2Bf^^^ for any value of r > 0. For A^^ = 0, the state of qubit B is not a separable state in the three qubit mixed 
state. 

For the state of Eq. (I26p . N^ {%) = for all the pure states 'pi appearing in the pure state decomposition. This is 
consistent with the earlier observation that the state does not have genuine tripartite entanglement. Genuine tripartite 
entanglement refers to tripartite entanglement due to correlations similar to those present in a three qubit GHZ-like 
state. Loss of a single qubit destroys this type of entanglement completely, leaving no residual entanglement. Only 
two qubit coherences are present in the state p^^^'^^[t) as evidenced by N^ {pi) = E2 {pi) for all pi. As such we may 
rewrite the Eq. (|49)) as 

NhDG{p)=Y.P^El{P^)■ (52) 

i 

Two qubit coherence can generate not only pairwise entanglement but also W-like tripartite entanglement. It is 
common practice to trace out subsystem A2 to obtain the entanglement of the pair AiB. One can, however, obtain 
a measure of 2— way coherences involving a given pair of subsystems {AiB or A2B or A1A2) from 2— way partial 
transpose constructed from the three qubit state operator p by restricting the transposed matrix elements to those 

for which the state of the third qubit does not change |12l |. For example, P2 ^~ ^^ is obtained from p by applying 
the condition 



(«1«2«3|P2 



" ^^ ■*''' b'l«2J3) = (il«2«3|p|«l«2J3); if ^ {^ - Si^Jm.) ^ '^, 



{■iii2i3\ P2^' ^'"^ \J1J2j3) = («i«2«3|pb'iJ2J3); for all otlicr matrix elements. (53) 

We also define the partial transpose P2^^"^^^^ 1 involving only the pair of subsystems A1A2, as 

3 
(iii2»3| P2^^'^'^^ \j1J2i3) = {jii2i3\ p\iiJ2i3) , if ^ i^ ~ '^i™j™) = 2, (54) 

m— 1 

and for all the matrix elements with 23 7^ J3 

{111213] P2^''^'^'' \hJ2j3) = {iii2i3\p\jii2J3) ■ (55) 

The partial transpose P2 ^ ^ is defined as 

3 
{iii2i3\ P2''"'''*' Vihh) = (iii2J3| P |iiJ2»3) , if X! (^ ~ ^in,-]^) = 2, (56) 

m— 1 

and for all the matrix elements with ii 7^ ji 

(ni2«3| P2''"''''' \i1J2J3) = (n«2«3| P IJ1J2J3) ■ (57) 
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FIG. 4: The negativity Npgj~,Q, Ep^j^ , and Ep^^p, versus t{= grjt) for s = 1.2. 



One can verify that 



and 



P2 =P2 +P2 -P' 



■^ A-\ -^ A-\ —BA-\ I -^ A-i ~ A-\ Ao -^ 

P2 = P2 + P2 - P- 



Using the Eqs. (|58p and ([55)1 along with the definition of £"2 in Eq. ([5^ we get 



^IsZ.g(p) = i?l5if^^(p)+£^f5D^^(p)> 



TV 



B 
P, 
A 
PSDG 



(P) 



^PSD yP)+^PSD \P)^ 



(58) 
(59) 



(60) 
(61) 



where partial negativity Epgj^ ^ {p) or Epgj^ ^ (jH) , associated with pure state decomposition, is defined as 

EfsE^ ip) -J2p^J2 ('^'J (p^)?-^^ K.) ' (62) 

Xfj^^i of the partial transpose {pi)^'^ ■ The partial 

\Tai-Aij 



and $ j is the /i*'' eigenvector corresponding to negative eigenvalue 
negativities -E'ps7,'^'^'(p) and Ep]fj^^^^(p) are given by 



E%lf^^^ {P) = E 



kH 



(Xp,«|(P»)2 



I -^/x,?/ ' 



where x~ ^ is the p*^ eigenvector corresponding to negative eigenvalue /3^.i of the partial transpose (^1)2^^ ■ 

The partial X— way negativities Epgp, ^ (p), Epgp, ^ (p) and -E2 ^^ ^ ^ calculated from pure state decomposition of 
state operator determine the pairwise entanglement of qubit pairs BAi, BA2, and A1A2. The entanglement available 



12 




FIG. 5: The negativity N^sdg, Epg^ and F{Wi) versus r(= grit) for s = 1.2. 



to Ai and i? if they have no knowledge of A2, as well as the probabilistic entanglement of subsystem BAi after a 
measurement has been made by A2, is determined by Ep^jj ^ (p). The eigenvalues and eigenvectors of p"^^"^^^ (t) are 
used to calculate, numerically, the partial negativities ^pgo^'"^^ (p), Epl;'^^'^ (p) and Epgj^^^ (p) = Epg^^^ (p). 



Figure (g]) displays iVp^ 



PSDG (^) 
B 



N 



PSDG 



(?)' ^PSD 



pAi-AiA2 
^PSD 

~r-i ^ ^ {p), and Epgp, ^ (p) versus parameter r for s = 1.2 and 
TT. The negativity N0gp,Q (p), and Epgj^ ^ (p) = Ep^p^ ^ (p) are plotted in figure ^ along with the fidelity 
F{Wi) defined as in Eq. ([33]). The maxima in F{Wi) plot are located in the regions with large Epgp, ^ (p) and finite 
Epg~p, ^ ^ (p), consistent with the fact that the reduced two qubit states p^'^i(r), p^^^{t), andp'^i^^(T) are entangled 
states. The partial negativities Ep^p, ^ (p), Ep^p, ^ (p) and Epg'p ^ ^ (p) detect the pairwise entanglement present 
in the composite state. Simultaneously nonzero positive valued Epg~p ^ ^, Ep]fp ^^ , and Ep^^ ^ (p), signal the 
W-like tripartite entanglement of three qubit state p'^^^^^{t). 

The calculated negativities Npgp^^ (p), Npgp^Q (p), Epg^ ^^ (p) and Epg^ ^ ^ (p) versus squeeze parameter s at 
T = 14.5 for 9 — TT, are shown in Fig. ©. The decrease in the probability of finding the cavities in vacuum state with 
increase in the value of s results in an increase in the probability of finding the system in an entangled state. The 
peak value of NpgpQ (p) occurs for s = 0.95 while NpgpQ (p) is maximized at s = 1.04. We have selected s — 1.2 
as the squeeze parameter for the reason that the probability of finding the qubits Ai, and A2 in entangled state is 
appreciable for this choice. 



VII. CONCLUSIONS 



We have studied the entanglement generation through interaction of squeezed light, shared by two remote cavities, 
with two atoms in cavity ci and a single atom in cavity 02- No direct interaction amongst the atoms takes place. 
Analytical expressions for the matrix elements of the three atom mixed state as a function of atom-field coupling 
strength and interaction time, are obtained. The decrease in the probability of finding the cavities in vacuum state 
with increase in the value of squeeze parameter s results in an increase in the probability of finding the system in an 
entangled state. Numerical calculations to study the dynamics of three atom entanglement show the qubit pair A1A2 
in Bell- like state entangled to remote qubit B. From the reduced state operators for three atom system, the linear 
entropy is calculated and compared with negativities of partially transposed state operators as well as negativities 
calculated from the pure state decomposition of the state operator. The bound entanglement of the mixed state is not 
detected by the global negativity, however, the negativity Npgp^Q (p) calculated from the pure state decomposition of 
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N' 



PSDG 



PSDG 

*- PSD 

cA,-A,A2 

^ PSD 



X =14.5 




FIG. 6: The pure state decomposition negativities A^pg^ig, Npg^Q, Epgj^ , and Epgj-, versus s for r = 14.5(= gr;t) 

the state operator detects the bound entanglement. The degree of entanglement depends on the squeeze parameter 
as well as the coefficient of reflection of the beam splitter. For r — 14.5, the pure state decomposition negativity is 
seen to be maximum for qubit Ai when s = 1.04, whereas the peak value of NpgjjQ (p) occurs for s = 0.95. A higher 
value of squeeze parameter, s — 1.2, has been selected to study the entanglement dynamics of the system since the 
probability of finding the qubits Ai and A2 in entangled state is higher than that for s = 1.04. The partial negativities 
calculated by selective partial transposition of the three atom mixed state detect the pairwise entanglement of qubit 
pairs AiB, A2B, and A1A2. The entanglement of three atoms is found to be W-like, no three qubit GHZ like quantum 
correlations being generated. It is well known, however, that a GHZ state can be distilled from a W-like tripartite 
state, by stochastic local operations and classical communication. Trapped atoms with long lived metastable states 
are favored as memory qubits. The advantage of having qubits Ai and A2 in Bell-like entangled state is that the 
entanglement between qubits Ai and B can be enhanced through local operations on qubit A2. The entangled states, 
generated in the scheme proposed here, can be used for implementing quantum communication protocols and for 
quantum computation. 
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